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Abstract 

The AGT conjecture relates M — 2Ad SUSY gauge theories to 2d CFTs. Matrix model techniques can be used to 
investigate both sides of this relation. The large N limit refers here to the size of Young tableaux in the expression 
of the gauge theory partition function. It corresponds to the vanishing of J7-background equivariant deformation 
parameters, and should not be confused with the t'Hooft expansion at large number of colors. In the first part 
of this paper, a saddle point approach is employed to study the Nekrasov-Shatashvili limit of the gauge theory, 
leading to define /3-deformed, or quantized, Seiberg-Witten curve and differential form. In a second part, this 
formalism is compared to the large N limit of the Dijkgraaf-Vafa /3-ensemble. A transformation law relating the 
wave functions appearing at both sides of the conjecture is proposed. It implies a transformation of the Seiberg- 
Witten 1-form in agreement with the definition proposed earlier As a side result, a remarkable property of M = 2 
theories emerged: the instanton contribution to the partition function can be determined from the perturbative 
term analysis. 
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1 Introduction 



A very powerful technique to evaluate matrix model integrals consists in expending correlators and partition func- 
tion in the matrix size N. Solving recursively Schwinger-Dyson, or loop, equations, the Eynard-Orantin topological 
recursion provides the expression of these quantities at all orders in the l/N"^ expansion [1]. The formalism survives 
the /3-deformation [2, 3] and can be applied to the Dijkgraaf-Vafa (DV) /?-ensemble representing Liouville correla- 
tors [4]. It was recently realized that this large technique can be employed directly within Liouville theory to solve 
perturbatively a set of Ward identities, leading to the semi-classical expansion of correlation functions [5]. In the 
spirit of AGT conjecture [6], a similar procedure should be available to treat directly the instanton partition function 
of four dimensional gauge theories with M = 2 supersymmetry.^ Such a technique would find appUcations in a 
large area of problems where localization is employed and results in a set of coupled integrals. 

The first step in this direction was made ten years prior to AGT by Nekrasov and Okounkov [21]. In this 
work, they studied the partition function of A/" = 2 (matter) SU (Nc) super Yang-Mills theories defined on the 
Q-background. They argued that when both equivariant deformation parameters ei, £2 tend to zero, the sum over 
Young tableaux involved in the expression of the partition function [22] is dominated by a partition of diverging 
size k ~ l/(eie2). This partition is characterized by a continuous density of boxes, the limit shape. Introducing 
a collective field action for the limit shape, they obtained its expression as a solution of a saddle point equation. 
They reproduced the Seiberg-Witten (SW) curve [23], and showed that the instanton free energy reduces to the SW 
prepotential. The strength of this technique lies in the possibilities of generalization. It has been used in [24] to solve 
the full set of quiver A/^ = 2 theories defined over a general (affine) Lie algebra [25]. 

Since the theory on is recovered as ei , €2 ^ 0, this limit will be referred as the SW limit. Following the 
AGT-DV correspondence, it is equivalent to the large (or planar) limit of the /3-ensemble at /3 = 1. Tuning /3 
with N, it is possible to define a double-scaled planar limit for the /3-ensemble which corresponds to the Nekrasov- 
ShatashviU (NS) limit of the O-background, namely €2 tends to zero while ei is fixed [26]. As in the SW limit, large 
N techniques can be developed to handle the partition function of M = 2 SUSY gauge theories [27]. Analysis of 
quiver theories was recently performed along these lines in [28], extending the work [24] to the NS background. 
One of the motivations for the present paper is to provide more details on this construction. To do so, we consider a 
slightly different point of view from [27], and take ei to be infinitesimal instead of finite. The main reason for this 
choice is to establish /3-deformed SW relations. 

The two celebrated formulas derived by Seiberg and Witten in [23] relate Coulomb branch vevs and prepotential 
to the cycle integrals of a differential form over an algebraic curve. Similar relations hold for the matrix model 
filling fractions and planar free energy, obtained as cycle integrals of the resolvent over the spectral curve [29]. 
The /3-deformation of these relations is very natural given the correspondence between SW theory and integrable 
systems [30]. This correspondence was originally observed for gauge theories defined on M'', and classical integrable 
models. It was later extended by Nekrasov and Shatashvili to quantum integrable systems and gauge theories defined 
on the J7-background with a vanishing parameter [26]. The non-zero parameter ei plays the role of a Planck 
constant in the quantized model. We focus here on gauge theories with SU{Nc) gauge group and Nj fundamental 
hypermultiplets. The associated integrable model is an Heisenberg XXX spin chain with Nc sites [31]. Remarkably, 
another integrable system shows up in the NS limit on the CFT side of the AGT correspondence. In the case of 
Liouville theory (Nc = 2), it is a reduced Gaudin model defined on a punctured sphere and characterized by a 
Schrodinger equation. The induced correspondence between quantum integrable models is referred as the bispectral 

'The original conjecture deals with gauge groups being a product of SU{2) groups and correlators of Liouville theory. Various verifica- 
tions of the conjecture can be found in [7, 8], and a proof was given in the series of papers [9]. The extension to SU{Nc) gauge groups and 
Toda field theory was proposed and investigated in [10, 11]. According to Dijkgraaf and Vafa [4], the conjecture can be reformulated using a 
/3-deformed matrix model. This proposal is discussed in [12, 13, 14], and also extends to higher genus and various gauge groups [15, 16]. The 
DV formulation of the conjecture has been demonstrated only for the case P = 1 which corresponds to the Hermitian matrix model [17, 18]. 

^The link between the moduU space of = 2 SYM instantons on R'' and Liouville theory was previously mentionned in [19], hinted by 
an earlier matrix model proposal [20]. 
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duality [32, 33]. The relation between Liouville theory and Gaudin model motivated the introduction in [34] of 
/^-deformed (or quantized) SW relations for the DV /^-ensemble. It was indeed observed that ei -corrections to 
the prepotential can be computed as Bohr-Sommerfeld integrals coming from the semi-classical treatment of the 
Schrodinger problem. These integrals extend the standard SW theory to the NS background. The fact that the /3- 
ensemble free energy obeys such deformed SW relations was later estabUshed in [35] using the Dyson collective 
field representation. The use of a semi-classical approach explains the choice of keeping ei infinitesimal. At ei 
finite, the algebraic SW curve is no longer well-defined, and the notion of "quantum curve" must be used [3]. 

This paper is divided into two parts. In the first section, taking the NS limit of the gauge theory instanton 
partition function, we derive the collective field action for the limit shape. In order to produce an equation of motion 
in agreement with the AGT correspondence, an additional cut-off term has to be included in the effective action, and 
we provide its expression. The resulting equation of motion was obtained in [27], it takes the form of a difference 
equation analogue to the Baxter TQ relation. We then proceed to the derivation of SW relations, and suggest ei- 
corrections to the standard definition of the SW differential form. The expression of the action cut-off term, and the 
induced SW differential are our main results. 

In the process of taking the NS limit, we observe the decomposition of the instanton effective action into the 
difference of two terms. Both terms consists of the same effective action, but evaluated for different collective fields, 
or densities. The first term describes the full partition function, and the density is obtained using the saddle point 
method. The second term involves a known density and reproduces the perturbative contribution to the free energy. 
Since the same action enters in both terms, and given the simplicity of the definition for the perturbative density, we 
conclude that the perturbative contribution to the free energy also determines the instanton corrections. 

The similarities between the treatment of the gauge theory partition function presented here and the standard 
large A'' approach to matrix models and /3-ensembles makes the comparison between gauge theory and DV model 
much easier, leading to a very explicit expression of the AGT correspondence. This correspondence is worked out in 
the second part of the paper for the case of an SU{2) gauge group and four fundamental flavors (SCQCD). It takes 
the form of a quantum change of variable analogue to the one used by Gaiotto to rewrite the SW curve [25]. This 
change of variable transforms the Schrodinger equation of the Gaudin model into the difference equation studied 
in the first section. It further implies that wave functions solving these two equations are related through a Mellin 
transform. This proposal, and its consequence on the deformed SW differentials, are investigated at the first orders 
in e. 



2 Gauge theory partition function in the NS limit 

In this section, we consider an = 2 gauge theory with SU (Nc) gauge group, and Nj massive hypermultiplets 
in the fundamental representation. When the number of flavors happens to be twice the number of colors, the /3- 
function of the gauge coupling r is vanishing, and the theory is superconformal. The dependence in the coupling 
constant will be encoded in the parameter q = e^'^*'^. For Nf / 2A''c, the gauge coupling is renormalized by a 
UV cutoff which generates an effective scale that we also denote q.^ The theory is considered on the Coulomb 
branch, and the vev of the adjoint scalar fields are denoted ai where I runs over the color index. Using a locaUzation 
technique, the instanton contribution to the partition function is reduced to a sum over coupled contour integrals 
[22, 36], 

^inst.feei,e„a,M = l + f;Cf^)' / {[^^^ {[ D[^, - cf^j). (2.1) 
The function D couples the integrals, it is a simple rational fraction involving only e-parameters, 

Dit) = (2.2) 



T 2N 

It is related to the parameter A[noj used in [21] by g = A, 



[NO] 



2 



The potential Q is a ratio of gauge and matter polynomials, 



M(f) ^= 
Q(') = AiFTT^TMify Mit) = Il(t-f^r), A{t) = ll{t-a,), (2.3) 

^ ' ^ ' r=l 1=1 

where denotes the mass of the hypermultiplets. The expression (2.1) is manifestly invariant under the sign flip of 

the parameters ai, fir, ei and £2, together with q — > {—l)^fq^ 

Contour integrals in the expression (2.1) of the partition function can be evaluated as sums over residues [22]. 
Poles on the real Une, namely for 0/ — 4>j equals to ei or €2, are avoided by a shift of the integration contours 
equivalent to the replacement ei^2 — ^ £1,2 + ^0. A similar treatment is applied to the poles at (pj = ai, with the shift 
a/ — a; + iO. Contributing poles are interpreted as the instanton positions, they are indexed by the boxes of Nc 
Young tableaux Yi, and will be denoted (pi^ij, where I is the color index, and {i, j) e Yi parameterizes the position of 

the box in the Zth tableau. We further denote ^ the height of the ith column of the Young tableau Yi, ni its number 
of columns and |y/| the number of boxes. ExpUcitly, the poles are given by 

(plAJ = ai + {i- l)ei + (j - 1)62, (2.4) 
with the indices (i, j) ranging from i = 1- ■ -ni and j = 1- ■ ■ xf \ The instanton partition function (2.1) rewrites 



/ 6 -j- 6 \ 

) = E n ^(^^ -MUn Qi^i)^ (2.5) 

{Yi,-,YM,}I,J=-i. 1=1 ^ ^^^^ ^ 

where the prime in the product denotes the omission of the factors with (pj — cpj = ±ea- To simplify notations, in the 
following the prime over products and sums will be omitted, and the suppression of diverging factors or summands 
always implicitly understood. Indices in the products run over the total number k = \ Yi \ of boxes in the partition, 
and the g-expansion is an expansion over the size of partitions. 

2.1 Sum over large Young tableaux 

The Young tableaux dependence in the summands of (2.5) can be encoded in two ways. In the work of Nekrasov 
and Okounkov [21], a shape function fy is associated to the Young tableau 5^ = (Ai > A2 > • • • > A„) as^ 



(2.6) 



/y(i|ei, €2) = - 2 1^1 - " ""'^"^'^ E [l* " " " ^^^^1 - \t - {i - l)ei| 

i=l 

where we introduced the convenient notation for the shift operator e'^^*g{t) = g{t + a). The function fy describes 
the profile of the tableau Y rotated such that it occupies the angular sector [7r/4, 37r/4]. Following [21], we assume 
that the Coulomb branch vevs are real, ordered as oi < 02 < • • • < a;, and separated by a distance greater than 
the Young tableau size e2A['^ x eim such that the support of the functions fy^it — ai), centered at t = a/, do not 
intersect. Cases where these conditions are not satisfied may be obtained from analytical continuation. We then 
define the total shape function / as the sum of Young tableaux profiles, 

f{t) = T.fY,it-ai). (2.7) 
1=1 



'^A different convention on the sign of the masses Hr is used in [6] and [22, 21]. For consistency with the next section, we employ here 

the convention of AGT [6]. 

^Our notation differs from [21], — 2/here = /[no]. Note also that the summation can be extended to i = 1 • • • oo if we set Ai>n = 0. 
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In the SW limit, the main contribution to the partition function is given by the limit shape /o obtained as the solution 
of a saddle point equation. In the NS limit with ei small, the Umit shape acquires ei-corrections that can also be 
derived from a saddle point technique, as will be shown below. 

On the other hand, Poghossian et. al. considered in [27] the density of instantons 

oo 

^mst.(0) = V S{ct) - 0,). (2.8) 

Our approach to the problem with ei infinitesimal provides a bridge between the formalisms used in [21] and [27], 
and we will work indifferently with shape functions / or densities p. 

As 62 -> 0, the size of Young tableaux diverges, and the height of each column is sent to infinity such that 
e2Xf^ remains finite. The row index e2{j — 1) becomes a continuous variable in the interval [0, €2Af ^] whereas the 
column index i remains discrete. Thus, the instanton positions 0;,^ j are now continuous variables (f)i^i that belong to 
the intervals [tj'j, ti^i] where t^^ = ai + (i — l)ei and ti^i = + \f^e2- In this process, an integral substitutes the 
sum over the index j, and for any function -F(0), 

In particular, the density of poles becomes a difference of sign functions, 

Nc rii 

Pinst. i4>)^-^J2Yl [^ig" - ^m) - Sign (t - ii)] ■ (2-10) 

1=1 i=l 

At this stage ni remains finite, it can be seen as a cut-off sent to infinity at the end of the computation. This is the 
point of view taken in [27]. In the NS limit, the shape function / is related to the instanton density as 

1 ^= 

fit) = 4rt.(*) + (1 - e-^^^*)pinst.(i). fm-it) = -2^1*- ^'1- ^2.11) 

1=1 

Sending further ei — > 0, we get upon integration the relation /o = /pert. + eiPinst. obtained in [21] (equ 4.37). 

The support of the density pinst. is rather compUcated, and it is easier to work instead with the densities associated 
to the bounds ■ and j of the integration over j, 

ptoit) = E E '^(* - *m)' pt(t) = E E - *'-^)- (2.12) 

1=1 i=l 1=1 1=1 

These densities are normalized to n = J2i and exhibit a natural decomposition over the color index. In the NS 
limit the triple summation over the indices I and {i, j) G Yi becomes a single integration involving the difference of 
densities pt - pto, 

Nc ^ 

^2E E ^('^MJ)^ / dtF\t)ipt{t)- Ptoit)), (2.13) 
1=1 {i,j)eY, 

where F^{t) is a primitive of F{t). The constant of integration is irrelevant here since pt and pto have the same 
norm. Upon integration by parts, the primitive F^{t) can be replaced by the original function F{t), and the dif- 
ference Pto — Pt by /Oinst. since p{J^^^ = pto — pt- The drawback of this method is the emergence of cumbersome 
boundary terms. We will not follow this path here, and simply work with the densities pt and pto- 
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We are now ready to take the NS limit of the instanton contribution (2.5) to the partition function. As €2 — > 0, 
we introduce the effective action 

k 

-5i„st.[{</)/}] = log (5Qo(0/)) + I E ^('^^ - '^•^)' (2.14) 

/ 7,J=1 

such that the instanton partition function (2.5) is approximated as £2 -> by 

^inst.(9,ei,e2,ai,/x,)~ J] 63 ^^'^"'e^*^'-'-. (2.15) 

The factor e^" is absorbed in the measure as the sum over Young tableaux is transformed into a functional integral 
[37]. Kernel and potential of the action are simply derived from the original expression (2.5),^ 

Applying the procedure detailed above to take the NS limit, we obtain 

>5inst. = / dt{pt{t) - pto{t)) f ds{pt{s)-pto{s))G^^{t-s)+ f dt{pt{t) - pto{t)) f log(gQo(i)) 

(2.18) 

with G^^ a double primitive of G, d^G^^ = G. The function G^^ (t) is regular at t = itei, however singularities 
may appear as we derive the equation of motion. The prescription for the regularization of poles at t = s ± ei is to 
shift ei — )• €1 + iO, as it corresponds to the case (pi = (pj ±. ei in the original expression (2.1). By definition, the 
instanton contribution to the free energy is 

•^inst. = - lim eie2log^inst. (2.19) 

e2->-0 

In the saddle point approach, it is obtained as Tm^t ~ — ei»Smst. where <Si„st. denotes the evaluation of the action 
functional (2.18) for the critical density, i.e. the density that extremizes the action. 

It is also possible to define correlators with respect to the instanton contribution to the partition function. 

In the NS limit, these correlators are approximated by 



e2/f;F(</>,)\ J^dtF\t)ipt{t) - ptoit)), (2.21) 

\l=i I inst. 



*The term coming from log D{i) has been symmetrized using (here (pij = (pi — 4>j) 

^ logD{(^u) ^ -eie2 ^ = -€162 ^ = T E ^2.16) 

In the process, we neglected the constant —tiG(ei) ~ — 2e2/ei coming from the subtraction of poles at (pij = ±£2. This constant is 
negligible in the NS limit as it just shifts the free energy by 2e|. 
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with a primitive of F, and the average performed with the critical density measure. In particular, using localiza- 
tion, correlators involving trace of powers of the adjoint scalar field $ reduces to coupled one-dimensional integrals. 
They can be computed from the generating function [27] 



iVc 



(tr e^*) = e""' - (1 - e"'^)(l - e'"') ( ^ e'^' ) . (2.22) 



1=1 \/=i 



To avoid complicated expressions in the following sections, we drop the index ei whenever the NS hmit is 
considered, and simply write e = e\. 



2.2 Revealing the perturbative part 

In this subsection, we reorganize the densities in order to extract the perturbative part, allowing us to consider the 
full partition function in the subsequent sections. To simplify our considerations, we first take the case of pure SYM, 
i.e. Nf = and M{t) = 1. In the saddle point approach, the Young tableaux columns (or the variables ti^i) are 
seen as random objects achieving a statistical equiUbrium. This equilibrium is obtained from the minimization of 
the effective action (2.18). From this point of view, they play a role similar to the eigenvalues of a random matrix 
of size n x nJ This matrix decomposes over the color index into blocks of size n^. The size of the blocks may also 
fluctuate, it is usually fixed by imposing filling fraction conditions in matrix models. The effective action 5inst. only 
depends on the difference pt — pto- The density pt involves the random height of the Young tableaux columns. 
On the other hand, the density pto contains known parameters, the Coulomb branch vev ai, as well as the "random" 
numbers of columns n;. Using the property 

(1 - e-^i^OPto(i) = Yl [^(* - - - - ^1"')]' (2.23) 
1=1 

we split the density pto = Ppert. — Pm into a known part ppert. and the random density Pm that satisfy respectively, 

Nc Nc 

(1 - e-^i^*)Ppert.(t) = Y^^^- ^"'i (1 - ^'"""'')Pniit) = ^ 5(i " a; - e^m). (2.24) 

1=1 1=1 

There is an ambiguity in the definition of these densities which takes the form of an e-periodic function. But this is 
actually irrelevant here as the densities will appear in the action (2.18) only through the difference p{t) — p{t + e) 
which is well defined. Moreover, this is not pertinent to our WKB approach with e infinitesimal for which the 
densities ppert. and Pm are infinite series of distributions. The main issue is the treatment of diverging moments: they 
must be regularized using the principal value at infinity, and the infinite contributions from p^^x. and pm integrals 
cancel each others. A more detailed discussion on these densities and infinities is provided in appendix B.^ 

The density p = pt + Pm contains all the random variables. Formally, i.e. up to a/ -independent divergences, it 
is possible to choose the following expression for p^^ , 

Nc oo 

Pni{t) = Yl Yl 6{t-ai-{i-l)ei). (2.25) 

1=1 i=ni+l 

^The method presented here is close but different from the one exposed in [38] and based upon the transformation of sum over partitions 
with Plancherel measure into matrix integrals [39]. The main difference lies in the definition of the eigenvalues since ai + e2{Xf^ — i + 1) 
is used instead of ti^i in [38]. It results in a model with a similar potential, but a different quadratic term. 

*The decomposition of the instanton density (2.8) pinst. = p — Ppea. into a perturbative density Ppen. = — Ppen. and p' = — p has already 
been introduced in [37], but without a clear interpretation. In the notations of [37], they split p = pq + pQ„ where pq, —pQg and p should be 
identified respectively with p, ppen. and pmst. here. The proper identification of ppen. can be found in appendix B. 
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Then, the density p plays the role of pt for Young tableaux completed with an infinite number of columns of zero 
height, 

Nc oo 

'"W = E E - °' - - 1)^1 - (2.26) 

1=1 i=l 

with ■* = for i > m. The decomposition of p in pt and Pm is related to the artic circle phenomenon [39], and 
Pm describes the frozen eigenvalues below the artic circle n^. 

However, we take here a different choice for Pm in order to keep a finite norm. Our choice of pn, is obtained 
perturbatively in ei in appendix B, it differs from the previous choice only by irrelevant "boundary" terms for the 
moments. The norm of the density ppen. is also computed in appendix B, and the computation is similar for p^ . The 
norm of pt is simply equal to n, and we deduce 

/ p{t)dt = [ Pp^n.{t)dt = --y]ai + l-Nc. (2.27) 

JR JR £ ^ 

For a gauge group SU{Nc), the sum over the Coulomb branch vev ai is vanishing. Nevertheless, it is instructing 
to keep it expUcit here. The first moments of ppert. and p are diverging, but the divergent terms compensate. Such 
boundary terms are independent of a, and only involve the cut-off and e. Thus, they do not contribute to the SW 
relations and can be discarded. We take 

I tppen.it)dt = -—y^ai{ai - e), / tp{t)dt = — \ 2eie2k - y^ai{ai + 2€ni) ] (2.28) 
Jr 2e ^ 7m 2e \^ ^ J 

where we used the vanishing of the a; sum to simplify the second formula. 

We now replace the difference of densities pt — pto in the action (2.18) with p — Ppert.- The main result of this 
subsection comes from the observation that the cross-term between p and ppert. cancels the potential term involving 
Qo for p. To make it appear, we should first fix the integration constant and define 

G^i{t) = (e^^* - e"^^*) i(log t-1), (2.29) 

such that both G{t) and G^^{t) are even functions. To evaluate the cross-term, we use a change of variable to transfer 
the shift operator from the kernel to the density Ppert., and then exploit the property (2.24). The resulting expression 
writes 

„ Nc I- ,.t 

/ dtdsp{t)p^n,{s)G^\t-s) = Y^ dtp{t){l + e'^') \og{s-ai)ds. (2.30) 

J'B? JR Jai 

Cancellation occurs with the potential provided we make a good choice for the integration constants, namely 

/ d.tp{t) [ \ogQo{t) = -Y] f dtp{t){l + e'^') [\og{s-ai)ds. (2.31) 

JR J i^-^ JR Jai 

This choice was left arbitrary as long as the integration constants were the same for p and ppert. ■ A similar property 
holds if we replace p with p^g.^ in the two previous formulas, it can be used reverse the sign of the quadratic term in 
Ppert.. In the remaining expression of 5inst., the densities are decoupled, 

«5inst. = <5sym[p] -«5sYM[ppert.], with Ssym[p] = -7; [ dtds p{t) p{s)G^ \t - s) + \og q [ dt {t + e/2)p{t). 

(2.32) 

We also fixed the integration constant for the primitive proportional to log q in the definition of 5sym- Again, this 
choice is arbitrary provided that we make the same one for p and ppert. • 
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It remains to interpret the action SsYM[Ppen.]- Discarding irrelevant boundary terms and divergences, the linear 
term proportional to log q reproduces the tree-level contribution to the free energy, 




(^l ^ tree-level — n 




1=1 



(2.33) 



To interpret the quadratic part, we need to introduce the function 7^ defined in the appendix and such that 



eG^\t) = ie'^ - e-^^*)(l - e'^)je{t). 



(2.34) 



Then, the quadratic term in ppgrt. reproduces exactly the one-loop vector contribution to the free energy, 




dtdsppen.{t)ppen.{s)G^^ {t - s) = — ^ [ye{ak - ai) + 7^(0^ - ai + e)]. (2.35) 

^ k,l=l 
k<l 



We have shown that the action 5sYM[Ppert.] supplies the perturbative (tree-level and one-loop) contribution to the 
free energy. Since by definition 5inst. provides the instanton contribution, the action >Ssym[p] = ^inst. + ^sYMiPpert.] 
generates the full free energy as T^f=2 = — e<SsYM[p]- As p is decoupled from Ppert., it is equivalent to extremize <Sinst. 
or 5sYM [p] to obtain the density p. Thus, the action 5sym [p] fully characterizes the partition function Z//=2 in the NS 
hmit. It is remarkable here that the very same action 5sym describes both perturbative and full partition functions. 
Knowing the definition (2.24) for Ppert., the action Ssym can be read from the perturbative free energy Jpert.- This 
definition reduces to a sum of delta functions centered at t = a; for p'^^^ in the SW Umit. Given its simphcity, it is 
natural to suppose that the definition (2.24) has a universal meaning. We conclude that the perturbative regime also 
determines the instanton contribution. It is of primary interest to figure out whether this is a general property of this 
set of J\f = 2 theories, and if it extends beyond the NS Umit. Unfortunately, such questions are beyond the scope of 
the present paper. 

In the presence of fundamental matter fields, the mass polynomials M{t) is responsible for a potential term in 
the action 5inst. given in (2.18). Since no new cross-term appear, the decomposition of 5i„st. into an action depending 
solely on p and the same action with p replaced by ppert. still holds. However, this action also acquires a potential 
term 



For a good choice of the integration constant, this new linear term in Sj<tf=2[ppert.] generates the massive contribution 
to the perturbative free energy at one loop. 



Consequently, the full free energy, including instanton and perturbative terms, is again given by the action depending 
only on p as 7^=2 = -e<SAf=2[p]- 

In [21], a trick was employed to re-introduced the mass term from the pure SYM action using a shift of the shape 




(2.36) 




(2.37) 
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function by an appropriate mass term. This trick generalizes to the NS Umit using the density Pmass defined as^ 

(e^a« _ e-^a*)p„,,,(t) = -J2S{t- lir). (2.39) 

r=l 

The action Sj^=2[p\ can be written as the difference of the SYM action 5sym[p + Pmass] ~ >5sYM[Pmass]- In this 
expression, the cross term between p and /Cmass generates the massive potential for p in the same way ppert. generates 
the potential term depending on the Coulomb branch vevs. The Unear term, and possible divergences, are eliminated 
by subtracting 5sYM[Pmass]- At several occasions we shall restrict ourselves to pure SYM, and this trick can be 
employed to recover the massive case. 



2.3 Shape function and SW limit 

In the limit e — > 0, it is simpler to work with the shape function / instead of p in order to recover the results of 
Nekrasov and Okounkov [21]. To relate these two quantities, first observe that from the definitions (2.24) and (2.1 1) 
of the perturbative density ppert. and the function /pen., we have the relation 

fv.n.{t) = -(1 - e-^^^*)pper,.(i)- (2.40) 

Differentiating the first expression in (2.1 1), using the relation between pinst. and pt, pto, and replacing the difference 
of these two densities by p — ppert., we end-up with 

r(t) = -(l-e-^i^*Hi). (2.41) 

This formula can be used to replace the shape function by differences of the density p. But it can also be inverted, at 
least perturbatively in e, to replace the density by the shape function. Considering the SYM effective action »Ssym, 
the shift operators coming from the kernel G can be transposed to the density with a change of integration variable. 
Using the previous relation (2.41), we obtain the expression of the action in terms of the collective field /", 

SsYuif"] = f dtdsf"{t)f"{s){l + e^^*)7e(t -s) + ^\ogqf dt t^f"{t). (2.42) 

As mentioned previously, matter can be reintroduced by adding a mass term to the function /, and subtracting the 
pure mass action. The norm and first moment of the shape function follow from the primary definition (2.7) of / 
[21], 

f . 

/ f"it)dt = -Nc, / tfit) = - V a;. (2.43) 

These formulas are in agreement with the relation (2.41) provided we take into account the boundary terms (see 
appendix B). 

In the SW Umit, the role of the density is played by — /q (i) = ^Po{t) where the subscript denotes the first order 
in the e-expansion. The RHS of (2.43) is independent of e, and /q also obeys these properties. The action (Ssym[/"] 
given in (2.42) reduces to 

5sYM[y^'] c.^J^^ dtdsf;;{t)f;;{s){t - sf (log (i-^) - , (2.44) 



'To make contact with the SW Umit, we may define a function /mass associated to pmass through a relation similar to (2.41), i.e. 

^/ 

(1 + e^^')/l3s(i) = ^S{t- Mr), (2.38) 

r=l 

In the SW limit e — >■ 0, we recover the definition used in [21] (equ 7.3). 
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which is indeed the effective action used in [21]. The solution to the equation of motion was found in [21] and will 
be re-derived in the next subsection. It exhibits a finite support F made of Nc intervals Ti = [a^~,a;^] such that 
ai G T;. On each of these intervals, — /q is positive if /o is concave, which may not be true everywhere. However, 
the positivity of the density is not required in the context of SW theory [29]. 

We will show below that each term of the e-expansion / = Yl'^=o ^^fn share the same support V as the limit 
shape /q . Due to the formula (2.41), it is also true of the terms in the expansion of p' . Moreover, both p and /" 
have a natural decomposition over the color index, and each component is vanishing on all but one interval T;. We 
deduce the "filling fraction" conditions 

I p{t)dt = --{ai-e/2), [ f"{t)dt = -l, [ tf"{t) = -ai. (2.45) 
2.4 Definition of tlie resolvents 

In the study of the planar Umit for matrix models, it is usual to introduce a resolvent associated to the density 
of eigenvalues. This quantity is known to have a branch cut on the support of the density, along which it has a 
discontinuity of — 2i7r times the density. For a large class of potentials, the resolvent satisfies an algebraic equation 
that defines the spectral curve. In the context of ^-ensembles, this algebraic relation becomes a first order non- 
linear differential equation, the Riccati equation. A wave function is frequently introduced to rewrite this problem 
in the form of a Schrodinger equation. Following this approach, we introduce here the wave function ^ such that 
dlogip{z) is the resolvent associated to p, 

logip{z)= / p{t)log{z-t)dt. (2.46) 
Jr 

Since we can equivalentiy work with the shape function, we also define a wave function u!{z) such that d log u!{z) 
is the resolvent associated to /", 

loga;(z)= / f"{t)log{z-t)dt. (2.47) 

JR 

These two quantities coincide with the ones employed in [27]. As a consequence of (2.41), the wave functions -0 
and oj are related through 

= (2.48) 

Placing the branch cut of the logarithm on M~, the function log V'(-z) is discontinuous on the interval ] — oo, ], 
and values above and below the cut read 

log ^(i ±iO)= / p{s) log \t-s\±iTr / p{s)ds. (2.49) 
Jr Jt 

Deriving this relation with respect to t, we recover the standard Feynman prescription for the resolvent, 

dt log 4j{t ± iO) = / ^^ds T iTTpit), (2.50) 

JRt — S 

provided we identify the principal value with the derivative of the weakly singular integral with kernel log |t — s|. 
The function has the same branch cut as its logarithm. On the other hand, a simplification occurs for u which 
has a branch cut only on the support F of /". The reason is that the discontinuity of logo; is constant outside of the 
support F, since its derivative is the resolvent of a density supported on F. Constants on each side of an interval F; 
differ by 2i7r times the integral of the density over this interval, which gives —1 in the case of /" (see (2.45)). There 
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is no discontinuity of log uj on the interval [a^^ , oo [, and these constants must be equal to a multiple of 2m. Upon 
exponentiation, we conclude that uj has no discontinuity outside of the support T. 

As noticed in [27], the wave function w is related to correlators with insertion of the adjoint scalar field 
Indeed, in the NS limit the correlator (2.22) can be computed using the formula (2.21). Replacing the difference of 
densities by p — /9peit.> the perturbative contribution cancels the first term in (2.22), and we end up with 

(tre^*) = - / dte^*f"{t). (2.51) 

Expanding in z, we deduce the moments 

(tr 1) = iVe, (tr = V a; = 0, = (tr = - dtff"{t), 



(2.52) 



where Ug with s = 2 - ■ ■ Nc parameterize the gauge theory vacua. Similar relations were obtained in [20] for the 
matrix model describing M = 2 SYM on M^. By identification of the terms in its z-expansion, we relate uj to the 
amphtude 

uj{zy^ = cxp (tr log(z - <!>)) ~ (det(z - $)) . (2.53) 
In this formalism, the Matone relation [40] is easily derived from the action (2.42), 

2.5 Equation of motion 

The critical density p minimizing the action Sj^=2 solves the equation of motion^^ 

^i^-Tl^ = ^ / dsp{s)G\t -s) = log {qM{t)) . (2.56) 

In terms of the wave functions ip or uj defined above, the equation of motion reads, 

gM(t)S^^ = l ^ qM{t)u:{t)u(t + e) = 1. (2.57) 
-^{t + e) 

Those equations were previously estabUshed in [27], they are valid for finite e and only holds at the discrete set of 
points t G In the case of an infinitesimal e, the support of the density condenses into the union of intervals T. 

In this limit, we need to introduce a regularization of the singularities at coincident instantons (pi = It results in 
the presence of an additional cut-off term in the effective action. Another cut-off contribution may also come from 
the replacement of the sum over Young tableaux with a functional integral over the density p. Such cut-off terms 
were first observed on matrix models by Dyson [41]. More recently, they have been considered in the context of the 
Dijkgraaf-Vafa /3-ensemble in [35].^^ In this model, the terms coming from the regularization of kernel and measure 
cancel in the SW limit where P = 1, and we recover the standard Hermitian matrix model. When /3 7^ 1, i.e. in the 
NS Umit, they should be taken into account as they reproduce the derivative term of the loop equations. 



''Working with the shape function /, we should instead consider 



Sf"{t) 

''See also [42] for a discussion in the framework of the 2D Coulomb gas. 



(l_e^^')at^%^=0. (2.55) 
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Our strategy will be as follows. We first concentrate on the SW limit, and suppose that there are no cut-off terms. 
We will derive an equation of motion valid on the whole complex plane, and recover the SW curve which validates 
our assumption. We then turn to the NS Umit, and consider the equation of motion derived from AGT conjecture in 
the next section. We deduce the expression of the cut-off term to be included in the effective action. This expression 
is further motivated by rewriting the kernel part of the action (2.42). 

In the SW hmit, we expand the kernel G^{t) ~ 2e/i and focus on the equation of motion 

2e-f ^ds = log (gMo(t)) , teV. (2.58) 

A possible e-dependence of the mass polynomial may appear in the presence of hypermultiplets in the anti-fundamental, 
and we denoted Mq its SW limit. It is more convenient here to work with the density /q = — epo, and consider 

2 / /^'(t) log \t - s\ds + log {qMo{t)) =0, teV. (2.59) 

The wave function ljq obeys the property (2.49) with p replaced by /q . Thus, the integral in the previous equation is 
simply the sum of the values of log coq above and below the branch cut. Exponentiating the equation, we deduce 

qMo{t)oJo{t + i0)ujo{t - iO) = 1, teV. (2.60) 

It implies that the quantity 

Po{z) = + qMo{z)iOo{z) (2.61) 
has no discontinuity on F. The asymptotic at z ^ oo can be derived from the relation (2.43), 

1 / 1 \ 

logu{z)^ -Nclogz + -Y,ai + 0\^-^j => Po{z) ^ z^^ +qz^f-^\ (2.62) 

Since there is no reason to expect any other singularity for Po(z), it should be a polynomial of degree Nc when 
Nf < 2Nc. In the following equations, we assume that Po{z) is a monic polynomial, and consequently when 
Nf = 2Nc Po{z) should be replaced by (1 + q)Po{z). The relation (2.61) between ooq and z defines a double cover 
of a sphere with Nf punctures at z = jir. This curve is identified with the SW curve of the gauge theory on 
[21]. This result justifies the assumption that no cut-off term is present in the effective action for the SW limit. The 
function ojq is multivalued and should be defined on a hyperelliptic surface of genus Nc with singularities at the 
poles of Mq{z). The two solutions of the quadratic equation (2.61) give the value of ojq on the two sheets covering 
the punctured sphere. 

The physical, or principal, sheet is by definition the sheet where ojq satisfies the asymptotic expansion (2.62). The 
value of uq on this sheet is given by Uq. The 2Nc branch points af are such that 

]\{z - a+){z - a-) = Po(^)' - 4gMo(z), (2.64) 
1=1 

which is indeed the equation (4.50) of [21]. In the superconformal case, the LHS must be multiplied by (1 — q) and 
Pq replaced by (1 -I- q)Po. 
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To go beyond the SW limit, we assume that the quantity P appearing in 



qM{z)uj{z + e)uj{z) - P{z)uj{z + e) + 1 = (2.65) 

has no branch cut over F. This is motivated by the study at finite e performed in [27]. However, a minor difference 
appears here in the sign of q. This sign flip is necessary to reproduce the correct SW curve as e ^ 0. It is also in 
agreement with the AGT conjecture, as shown in the next section. The arguments employed on Pq can be repeated to 
show that P{z) must be a monic polynomial of degree Nf., and has to be replaced by (1 + q)P{z) when Nf = 2Nc. 
Trading w for the ratio (2.48) of ip functions, the equation (2.65) produces a difference equation analogue to the 
Baxter TQ relation, 

qM{z)^{z - e) - P{z)'^{z) + ^{z + e) = 0. (2.66) 

It now remains to relate the condition of vanishing discontinuity for P{z) to the equation of motions. Our claim 
is that this is implied by the equation of motion (2.56), provided we correct it to account for the cutoff term in the 
collective action. The requirement Disc P = on T is equivalent to the condition, 

q-'M^u+u~utu;;F = l with ^ ^ (^"^ " ^~)' (2.67) 

u!+ui~ (a;7 — LJe r 

where oj^ is a shortcut notation for the function with shifted argument uj{t + e), and the superscript it denotes the 
position with respect to the branch cut, uj^ = uj(t± iO). The principal value regularization of the equation of motion 
derived from the action 5a/=2 produces the same identity with P = 1. The multiplicative alteration of the equation 
of motion is the sign of an additive correction to the action. Since a cutoff term arises from the regularization of 
kemel and measure, it must be independent of q or M, which is indeed the case of the function P. In addition, 
cut-off terms usually involve the density by itself, and not the principal value of the resolvent. It implies that P must 
depend only on the ratio u!~^/oj~, possibly shifted by e, and cannot depend on the product w+w-. The expression 
we found for P again satisfies this property. Put it otherwise, the cut-off terms we add here cannot be replaced by a 
modification of kemel or potential. 

The precise form of the function P can be given in terms of the shape function /, 

nt)= 'fli'}^j , logP(i) = 2(l-e^^*)log|sin(7r/'(i))|. (2.68) 

sm^ Trj'{t + e) 

We notice that P is real. We also verify that as e — > 0, P = 0(e) gives no contribution in the SW Umit. In general, 
at the order 0(e") where the equation determines f", the cut-off term involves only shape functions of a smaller 
order The corrected saddle point equation reads 

(1 + e'^*) / f"{s) log |t - s| + (1 - e'^*) log | sin(7r/'(t))| + loggM(t) = 0. (2.69) 
Jk. 

To motivate the expression for the cut-off term, we may investigate the quadratic part of the action. Since the 
potential plays no role in this discussion, we choose here to consider the SYM action with M{t) = 1. To reproduce 
the amended equation of motion, the collective field action Ssym should be modified to include the term 

<5sym[/"] = -- / dtdsf"{t)f"{s)%{t -s)+ [ l\f{t))dt + 1 log g / ef"{t), (2.70) 

with the function l\x) such that l\x)' = log | sin(7rx)| and the kemel 7^, defined in appendix A, being a com- 
bination of 7e functions. The origin of the cut-off term can be understood by noticing that unhke the kemel G^^ 
associated to the density p, the kemel 7^ is not an even function, but satisfies 

%{-x) = %{x) - eH'ix/e). (2.71) 
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This allows to rewrite the quadratic term of the action as 



-- / dtdsf"{t)f"{s)%{t -s) = -- f dtdsf"{t)f"{s)%{\t -s\)-e f dtdsf"{t)f"{s)l\{t - s)/e). 

(2.72) 

The second term in the RHS is weakly divergent for t = s, its regularization produces the cut-off term in the 
action (2.70) with a cut-off length ~ ef{t). As for the /3-ensemble, a more rigorous derivation may be obtained by 
considering loop equations. However, those have yet to be established for the present model, and we postpone this 
discussion to a forthcoming publication. 

We conclude this discussion on the equation of motion with a few more comments. The equation (2.65) has to be 
understood as a series expansion in e for the quantities cj, P and M. At the first order, we recover the two solutions 
(2.63) for ujQ. At the next to leading order, it provides the correction 

qMiujl + gMowo^o - Pic^o - Pquj'q 
= 2,Moa;o-Po ' ^'-''^ 

where the subscripts denote the order in the e-expansion. This function oji shares the same branch cut as and 
its derivative. It is actually easy to see that at the order 0(e"), the equation is always Unear in ujn and involves no 
derivative of this quantity. Thus, is simply determined by the lower order functions 0Jk<n and their derivatives, 
together with and M^. We conclude by recursion that all functions a;„ have the same branch cut as wq. They 
may however involve branch points of higher order, but always of the square root type, i.e. (z — af)"'/'^ with n^'L. 
Consequently, all orders corrections /„ to the limit shape, and p'^ to the density p' after inverting (2.40), share the 
same support F. As for the /3-ensemble, this may not be true anymore when we consider the resummed series at 
finite e [3]. 



2.6 Deformed Seiberg-Witten relations 

We now proceed to the derivation of e-deformed SW relations. We follow the approach used by Chekhov and 
Mironov in [29] on the Hermitian matrix model. This derivation was later extended to the Dijkgraaf-Vafa ^-ensemble 
in [35]. The major difference is the presence of a cut-off term in the effective action close to the one obtained in 
(2.70). This term results in poles contribution at the branch points which modify the SW relations. From this point 
of view, the treatment presented here is very similar to the case of /3-ensembles. 

The dependence of the collective action in the Coulomb branch vev ai is encoded in the density, and the derivative 
of the action with respect to ai is simply 

^'^•^=^ - ^ dAt)C{t) where C{t) = (2.74) 



dai 7k dai dp{t) 

As a consequence of the equation of motion (2.56), C{t) is a constant q on each interval F;. In the context of 
(deformed) SW theory, these constants should be allowed to be different [29].^^ The remaining integrals can be 
evaluated by differentiating the relation (2.45), leading to 

= -ci/e. (2.76) 



da, 



'^In the SW limit, ttie constants c; are related to the Lagrange parameters ^; introduced in [21]. This can be seen after writing C{i) with 
the help of the limit shape fo{t). Taking here M{t) = 1 for simplicity, the condition 

C{t) = 2 ^ dsf(;{s){t - s) (log 4^ -i)=ci, for teVi, (2.75) 



V[NO] 



is the analogue of the equation (4.40) in [21]. 
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The constants q, or more precisely their differences, can be obtained as cycle integrals of a 1-form over the SW 
curve. To do so, we need to define a function C{z) on this curve, i.e. a bi-valued function on C with a branch cut on 
r. We require that the sum of the values above and below the branch cut F; is equal to 2ci. We also introduce the 
cycles Ai that circle the branch cuts Ti in the physical sheet, and the dual cycles Bi that satisfy Akr\Bi = Sk,i where 
indices range from 2 to Nf.. The cycles Bi relate the two sheets by passing through the first and Ith branch cuts. The 
difference of constants q can be expressed as a cycle integral over Bi^i — Bi which is passing through the Zth and 
(/ + l)th branch cuts, 

dS^=2 dSM=2 1 ^ ^^^^^^^^ Z = l...iV,-l, (2.77) 



dai dai+i 2e Jbi+^-Bi 

with the trivial cycle Bi. In this expression, dC can be shifted by any meromorphic function v{z) on C provided it 
has no poles inside the integration contour, i.e. at the branch points. This property, which also holds for Ai-cycle 
integrals, can be used to simplify the definition of the SW differential. We will set 2dS = dC + v{z)dz with v{z) 
chosen appropriately. Taking into account that for an SU {N) gauge group, the sum of a^-partial derivatives of the 
free energy is vanishing, we deduce the SW relation 

= / dS. (2.78) 

oai Jbi 

We now determine C{z). Since we are only interested in the differences q — it is simpler to directly look 
for dC{z). It is also more convenient to work here with the limit shape / instead of the density p. The expression 
of C{t) is obtained by evaluating the variation of the action with respect to 6p{t), and then replacing the remaining 
densities by shape functions using (2.41). In this derivation, we have to take into account the cut-off term given in 
(2.70), we get 



Hit 



(1 + e'^*) / /"(s)log |t - s| + (1 - e'^«) log | sin(7r/'(i))| + \ogqM{t). (2.79) 



It suggests to take 



dC{z) = (1 + e''') log uj{z) + (1 ^ e''')L{z) + log qM{z), (2.80) 
with a function L{z) defined on the SW curve, and satisfying for t G T the relation 

L{t + iO) + L{t - iQ) = 21og I sin(7r/'(i))|. (2.81) 

At the origin of the function L{z) is the presence of the cut-off term in the action. A similar function, also denoted 
L, appeared in the treatment of the DV ^-ensemble [35]. It can be expressed with the help of another function b{z), 
also defined on the SW curve, and with a vanishing real part on F, 

L(z) = log6(z), with 6(t ±iO) = ±zsin7r/'(t) for t G F. (2.82) 

The function h{z) is a solution to the quadratic equation 

,2 l(a;+-a;-)2 



(2.83) 



where the functions u}^{z) are the values of a; on the two sheets. They can be obtained recursively using the equation 
(2.65) and starting respectively with the functions uj'^iz) found in (2.63). We observe that the discontinuity of L 
over F; is a constant, and it gives no contribution to dC. However, to suppress the function L in the definition of the 
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SW differential dS, we further have to show that it has no pole at the branch points. It turns out that such poles exist 
at the first order, 

(1 - e^^)L(^) ~ ^51ogMo(^) -'-J2 (2-84) 

This term is responsible for shifts of — e/2 for the A; -and i?/ -cycle integrals. It is possible to show that there are 
no pole contributions at higher orders. To simpUfy later comparison in the context of AGT conjecture, we neglect 
here the function L in the definition of the SW differential, and take 

dS = ^ \og{oj{z)oj{z + e))dz. (2.86) 



The counterpart is the presence of the e-shifts in the relations. 



dai jj^^ 



= (f dS-me{l-l), l = 2---Nc. 
Jb, 



(2.87) 



Upon integration by parts we may write 



dS = ~zd\og{u;{z)u;{z + e)) (2.88) 



and we recover in the SW limit the usual 1-form dS = —zdu)/u. 

Once the differential form has been determined, it remains to verify the A/-cycle relation. Since 5 log a; is the 
resolvent associated to the density f"{t), and given the properties (2.45), we find 

^^^dS = -^J^ + e-'^')tdt = ai-^. (2.89) 

We recover here the e correction that was first observed in [43] in the context of the Dijkgraaf-Vafa matrix model. 
It is the consequence of poles at the branch points, and can be absorbed by including the function L within the 
definition of the SW differential. 

In the next section, we investigate a transformation of the SW differential into the equivalent 1-form for the DV 
^-ensemble. It will be simpler there to consider the wave function ^ instead of oj, 

dS = ^z{e'^ - e-'^)d\ogi>iz). (2.90) 
We claim that upon a shift of integration variables this differential is equivalent to 

dS = -ed\og'ip{z)dz. (2.91) 
It is easy to check using (2.45) that it gives the same result for the ^^-integral. 



'^The proof goes as follows. We denote T[R] the ring of functions of the form A + BR where A and B are meromorphic and R is the 
square root of the polynomial given in (2.64), with branch points at q^. This set is stable under multiplication, ratio, and differentiation. 
We easily check that uio, and by recursion all uin belong to this set. It is also true of 6^ at any order. Expanding the logarithm, we find that 
L = log bo + J^[R] where 6o is the first order in the e-expansion of the function b, and J^[R] indicates the presence of a function in this set. 
We deduce 

(l-e'^)L(0) = -edlogboiz) + dT[R]. (2.85) 

The first term generates the 0(e) correction presented in (2.84). The second term designates the derivative of a function in J^[R]. It has no 
residues at the branch points and gives no contribution to cycle integrals. 
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3 AGT relation in the NS limit 



In order to compare our results with the DV /^-ensemble [4], we restrict ourselves to the gauge group SU (2). Our 
aim is to relate the density of eigenvalues of this model, or more precisely the associated wave function to 
the corresponding quantity on the gauge theory side. Higher rank gauge groups can be considered at the price of 
introducing a /3-deformed matrix chain [15]. These models involve several eigenvalue densities. Since only a single 
density is needed for the Nekrasov instanton partition function, the matching becomes more involved. On the other 
hand, the quantum change of variable presented here easily generalizes to such more complicated cases. 

For simplicity, we will only treat the case of four flavors: two in the fundamental representation of the gauge 
group with masses n±, and two in the anti-fundamental with masses jl±. Models with a smaller number of flavors can 
be obtained by taking a suitable limit [44]. Since the equality Nf = 2Nc is satisfied, the model is superconformal. 
It is sometimes referred as superconformal QCD, or SCQCD. The gauge and mass polynomials entering into the 
expression (2.1) of the instanton contribution to the partition function are respectively of degree four and two, 

M{t) = {t- n+){t- n-){t- e + fi+){t- e + jl-), A{t) = t'^ - a'^ . (3.1) 
3.1 Dijkgraaf-Vafa /3-ensemble and Liouville theory 

We start from the Liouville correlator of four vertex operators inserted at the points 0,1, q and oo. Integrating over 
the Liouville zero mode, the correlator exhibits a pole whenever the neutrality condition upon charges of the vertex 
operators is satisfied [45]. The residue can be written in the form of Selberg integrals similar to those appearing in 
the work of Dotseenko and Fateev [46]. Upon the non-trivial assumption, investigated in [47], that the prescription 
of the integration contours for screening charges is equivalent to the filling fraction condition in the standard large 
N approach to matrix models, we are led to study the ^-ensemble partition function 

Zp{a,mi,g,/3)= / He-^'^^'UXi [J {\i - Xjf , (3.2) 

1=1 i,J=i 

with g'^ = —£162/4 and /? = —€2/(1- The Penner type potential is reminiscent of the vertex operators. It introduces 
singularities at x = 0, 1, g with residues being the rescaled Liouville charges mo, m\,m2, 

= (mo + ^)logx + milog(x- 1) -|-m2log(x-g), ^=^(ei + e2). (3.3) 

We choose here the potential corresponding to the model studied in [48, 43, 35] and mostly keep the notations of the 
reference [35]. The dependence on the fourth charge moo appears through the neutrality condition 

mf + 2^gN = 0, (3.4) 

/=0,l,2,oo 

where N is the number of eigenvalues. According to the AGT correspondence [6], Liouville charges are related to 
the mass of hypermultiplets through 

A*± = w-2 ± mQ, /i± = mi ± moo- (3.5) 

The NS limit is equivalent to a semi-classical limit for Liouville theory with heavy vertex operators. In this limit, 
the Liouville charges, as well as the background charge, have to be rescaled by a factor ~ g, and the bare charges 
tend to infinity as 5 — )• 0. The number of eigenvalues N is also sent to infinity while g^/P ~ 62 goes to zero, such 
that the neutrality condition (3.4) is satisfied. In this limit, the planar free energy of the ^-ensemble defined as 

= limg'^logZp (3.6) 
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is related to gauge theory free energy Jvv=2 through [43] 

Tm=2 = 4J> + ((mo - m2f + 4h{mo - K)) log q - 2(mi - 2h){m2 - 2h) \og{q - 1) + est, (3.7) 

with h = e/2 and a q- and a-independent constant. As in the previous section, a density can be introduced, now 
associated to the eigenvalues A/, 



limg,/p(j25i\-Xi)y (3.8) 



P/3(A) 

This density extremizes the Dyson effective action [41] 

Ap] = I V{X)piX)dX + J log |A - fi\p{X)p{ti)dXdii + hj p{X) log p{X)dX, (3.9) 

where the last term in the RHS is a cut-off term coming from the measure and the regularization of the kernel at 
coincident eigenvalues. The equation of motion derived from this action is equivalent to a Riccati equation for the 
shifted resolvent H{x), 

Hixf + ed^Hix) = T(x), Hix) = 2 [ + v'(x). (3.10) 

J X - X 

Apart from masses m/ and equivariant deformation parameter e, the function T{x) also depends on a variable E. 
This variable is roughly equal to the gauge scalar amplitude U2 = (tr <^>^), up to a gauge coupling dependent factor 
and a translation involving masses [49]. This amplitude is defined in (2.52) and parameterizes the gauge theory 
vacua. It can be evaluated in terms of the Coulomb branch vev a using the ^-cycle SW relation, equivalent to the 
filling fraction condition for the /3-ensemble. With a slight abuse of terminology, the function T{x) is referred to 
as the "stress energy tensor" since it plays the role of this quantity in the Coulomb gas description [50]. It presents 
double poles on the sphere CPi at x = 0, 1, oo and reads 

T(x) = + (3.11) 

x^{x — 1)^(2: — qy x{x — l){x — q) 

where P4{x) is a polynomial of degree four in x, quadratic in e and independent of E. Its expression will not be 
given here, but can be found in [35] (equ. 2.20). In the SW hmit, the Riccati equation (3.10) reduces to an algebraic 

relation between H and x that defines the spectral curve. This spectral curve is a double cover of the sphere with four 
punctures, it has been identified with the SW curve in [48, 49] . Introducing the wave function ^(x) as H = ed log ip, 
the Riccati equation is equivalent to a Schrodinger equation, 

{y^-T{x))^jJ{x) = 0, y = ed^, [y,x] = €. (3.12) 

Here, a notion of quantum curve emerges as the operator acting on ip is simply the spectral curve of the model at 
h = 0, with the function H replaced by the operator y [3]. This Schrodinger equation can be derived directly within 
Liouville theory, without the ^-ensemble as an intermediate step [51]. 

SW relations for the DV /3-ensemble were derived in [35] using the Dyson effective action. It was argued that 
because of poles at the branch points, it is necessary to subtract to the usual SW differential dS = H{x)dx a function 
related to the cut off term of the action. It was then defined as dS = dxTjdx with dxT] solution of 

H{x) = dxtiix) - ^dxlogdxr]{x). (3.13) 

However, to simplify the comparison with the gauge theory, we keep here the definition dS = H{x)dx of the SW 
differential, and instead consider the e-corrected relations 

— I dS = a- e/2, <l> dS = 4^ + me. (3.14) 

2z7r /_4 /e da 

These formulas are compatible with those found in (2.89) and (2.87) for JV=2 under the ACT relation (3.7). 
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3.2 A quantum change of variable 

We now express the AGT relation in the NS Umit using our formalism. Similar approaches can be found in [32, 33], 
where the connections with integrable models is discussed in more details. In particular, the quantum change of 
variable we present here has been found in [33]. We reproduce it for the sake of completeness, and present a slightly 
different derivation. To keep expressions simple, we expose in this subsection only the massless case Ai± = A± = 
0. The massive case is treated in appendix C and requires extra care to cancel poles. Our starting point is the 
Schrodinger equation (3.12), 

2^(6/2)2 {q-l)E 



x{x — l)(x — q) ^/ ^ ^ 

In order to emphasize the fact that we are performing a quantum computation, we used a bracket notation, and 
commutation of variables should be performed cautiously. The first step is to multiply on the left by the product 
x{x — l){x — q) in order to cancel the poles a.tx = l,q and reduce the order of the pole at x = 0, 

{x\{x - - ^)x-^(x\2 + (e/2)2) -{q- 1)E\'iIj) = 0. (3.16) 

Then, we introduce the quantum change of variable z = xy — e/2. Commutation relations between x and z are 
derived from the canonical commutator between x and y, 



\x, z] = —ex =^ xqiz) = qiz ~ e)x, 

_ ^ I ^ ^^^^^ 



~^,z\ = ex ^ ^ X ^ g{z) = g{z + e)x ^. 



They imply that x acts as a shift operator on functions g{z). These commutation relations remain unchanged when z 
is shifted by an arbitrary function of x. This property is exploited in the massive case to further define z = xy + s(x) 
with an appropriate function s(x). Expressing y in terms of z in the previous equation (3.16), we get 

<x|(x - l)(x - q)x-^z^ - (g - l)E\il)) = 0. (3.18) 

Classically, this equation is simply the Seiberg-Witten curve [23] in the Gaiotto coordinates [25]. We now define a 
representation of the operator x in the basis \z > compatible with the commutation relations (3. 17). This requirement 
ensures that the commutator between x and y is preserved so that the change of variable (x, y) {z,p) withp = edz 
is canonical. We choose here 

X = r{z)e-P, x-^ = ePr{z)-^ = r{z + e)-^eP, with r{z) = qz^ . (3.19) 

The operator x acts in the z basis as a shift operator with an additional multiplicative factor r. It transforms the 
Schrodinger equation into a difference equation, 

-edz ^^ , „^D/.^ , „e9zl„/.\_n 2/ ,^2 r,! ^\ 2 l"? 



x|gM(z)e-'^^ - (1 + q)P{z) + e^''^ = 0, M{z) = z\z - ey, P{z) = z^ - J^p^^- (3-20) 

which is precisely the relation (2.66) studied on the gauge theory side. We also recover the proper expression for the 

mass polynomial M{z) as given in (3.1), and derived the expression of P{z) which is indeed a monic polynomial 
of degree two. In the massive case, the functions s(x) and r{z) involved in the change of variable acquire a mass 
dependence, 

s{x) = n+- xV'{x), r{z) = q{z- n+){z- nJ), (3.21) 

where the /3-ensemble potential ^(x) is defined in (3.3). This change of variable is not the only one possible. 
Other choices include: sign flip of y in the definition of z, translation of s{x) by a constant, flip of the masses 
mi 2 — )■ e — mi 2 appearing in the residues of s(x) at x = 1 and x = q,... These different choices lead to different 
mass polynomials. However, it can be shown that they describe the same theory by exploiting various invariance 
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such as the sign flip of a, or the exchange of fundamental and anti-fundamental hypermultiplets. Although probably 
worth of investigation, the inventory of all possible choices will not be done here, and we focus on the one given in 
(3.21). 

We conclude that the wave functions ^(x) and ^(^{z) represent the same state \ip > expressed either in the \x > 
or \z > canonical bases. It means that the two wave functions are related through an integral transform with kernels 
k{z, x) and k~^{x, z). Introducing the representation of the operator z and x in the \x > and \z > basis, we find 
two constraints on the kernels which are solved by 

A;(z,x) = a(2)x-^/^ei^''^(^)'^^°sa;, k'^x, z) = k{z,x)-\ and -^^^ = r{z). (3.22) 

a{z - e) 

Like the function lip, these kernels are determined up to an e-periodic z-factor. This degree of freedom arises on the 
gauge theory side because the mapping {y, x) — )• {p, z) is not one to one, but leave the possibility to shift ^? by a 
multiple of 2z7r. We also verify that formal inversion properties hold, 

/dlogx oc - .'), /d. c i(log. - log.'). (3.23) 

Thus, we have shown that the AGT relation takes the form of a MelUn transform between wave functions, 

ip{z) = a{z) / dxx— ^e--^^''>ip{x), V(^) = / dzx^ aizy^ipiz). (3.24) 

Jo 2z7re 7c-ieoo 

The appearance of the MelUn transform was expected here since when r{z) = 1 the momentum p is simply minus 
the logarithm of the coordinate x, and the change of canonical basis between p and z is known to be characterized 
by an exponential kernel. It is worth noticing that in (3.24) the potential term compensates the translation by V' 
in the definition (3.10) of the resolvent H. This shift also appears when the wave function ip is expressed as the 
^-ensemble correlator 



1p{x) = 62^ 



lim^nix-A,)^ (3.25) 



The function a{z) plays the role of a potential term but for the function -0(2;), it can be taken to be 

a{z) = (ge2)-/T [1 + {z - /i+)/e] T [1 + (z - . (3.26) 

Replacing — (T~^^ boils down to multiply u hy r which corresponds to shift the density f"{t) by mass delta 
functions, as can be seen comparing (2.47) with 

log riz) = V / ,5(i - /x±) log(2; - t)dt + log q. (3.27) 

_l_ JR 

As for the \x > space, there is a choice of coordinate, through r{z), such that the difference equation over co becomes 
trivial. However, finding this coordinate is of the same difficulty as solving the full equation. The convergence of the 
first integral transform in (3.24) depends on the behavior of the function ipi^) at zero and infinity. The Schrodinger 
equation (3.12) has two linearly independent solutions 'tp± with a different asymptotic, 

^±{x)e-eVi^) ^ x-e(^'++^^±\ x ^ oo. (3.28) 

The physical solution, or more precisely the wave function reproducing the resolvent on the physical sheet, is 
As z — > 0, the two solutions behave as 

V'+(x)e-^^(^) - 1, V-(a;)e-^^(^) x^O. (3.29) 
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Supposing Re e > 0, the integral over z converges when n± > z > —jj,± for respectively the wave functions 'ip±{x). 
Values of ip{z) for z outside this interval must be obtained from analytical continuation, which implies a deformation 
of the integration contour in (3.24). The transformation is worked out in details in appendix D for the simple case of 
a stress-energy tensor with only one double pole. 

In the limit e — 0, the inverse Mellin transform can be approximated by the method of steepest descent. In 
general, employing this saddle point technique on a Laplace transform leads to a Legendre transform at the first 
order of approximation [52]. For this purpose, we introduce a notation for the phase of wave functions, 

= ei-^^^), = e--J^'\ (3.30) 

Using the CFT methods for matrix models [50], with a natural generalization to ^-ensembles by turning on a back- 
ground charge [15], the function (/)(.t) is identified with a Coulomb gas field. This field is roughly the chiral part of 
the Liouville field after integration over the zero mode and a Wick rotation [35]. As e ^ 0, the transformation 



, , c+ieoo 



"I /*C-pl.COO 

= / (3.31) 

2i7re Jc—ieoo 



reduces to a Legendre transform between 4> and cf), 

(p{x) — V{x) = (z — log X — {4>{z) + e log (t{z)). (3.32) 
The mapping x{z) is obtained as a solution of the saddle point equation 

logx = — e51og(-0(z)/o-(z)). (3.33) 

In the SW limit, this is equivalent to set x = r{z)ujQ{z). Up to the function r which can be absorbed in a redefinition 
of cj, we recover the classical change of variable performed by Gaiotto in [25] : the original differential form zdooo / ujq 
is replaced by ydx with x = ooq and z = ujoy = xy. Thus, the transformation {p, z) — )• {y, x) can be interpreted as 
the quantized form of the Gaiotto change of variables. 

On the Liouville side, the SW differential is simply dS = dcf). It is actually more convenient to write it as 
dS = xdx4>{x)dlogx. The Legendre transform imphes 

diogxHx) = xV'{x) + {z- fx+), (3.34) 

and log a; is given by (3.33). The terms containing V', /Lt+ or a generate no contribution to A- or S-cycIe integrals 
and can be neglected in the definition of the SW differential. They will however modify the pole structure of the 
1-form, replacing the poles at x = 0, 1, (7, oo by poles at 2; = /x±, /i±. We deduce the expression of dS in the basis 
\z > at the first order, upon integration by parts 

dS = ed log tp{z)dz = -d^, (3.35) 

i.e. dS = zdlogojo. We recover here the expression (2.91) found in the previous section, up to a sign. This sign 
seems to indicate that the physical and non-physical sheets are interchanged under the transformation (3.24): the 
physical wave function ijj^ gives rise to which produces the values of the resolvent on the non-physical sheet 
(a;+) through the relation (2.48). 

At the next order, we need to take into account the square root of the Hessian in the evaluation of the integral, 
and the Legendre transform (3.32) get e-corrections in the form, 

(^(x) - V{x) = iz- ix+) log X - {4>{z) + e log a{z)) - (e/2) \og{d^ log(^(z)/a(z))) (3.36) 
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However, the saddle point equation (3.33) and the relation (3.34) still holds. It implies that the SW differential is 
still given by the expression (3.35), in agreement with the results from the gauge theory. Higher order corrections 
to the saddle point approximation do not modify the equation (3.34), and the result (3.35) actually hold at all order 
in e. It can be seen as a proof of the matrix model version of the AGT conjecture in the NS limit and up to an 
a-independent term. To go further and prove the AGT relation (3.7) including the o-independent constant requires 
a better understanding of the transformation for the densities and actions. Finding the correct change of variable 
for the functional integral of the collective field theories would be a first step in the derivation of a matrix model 
transformation between the /3-ensemble and the gauge theory partition function. 

4 Summary and discussion 

In this paper, we employed a saddle point technique to study the partition function of = 2 4d gauge theories in the 
NS limit. Starting from the instanton partition function, we introduced the densities referring to perturbative and full 
partition functions. The later one can be treated as a collective field minimizing an effective action. To reproduce 
correctly the results induced by the AGT conjecture, it is necessary to include a cut-off term in this action. This 
term is responsible for an additional contribution to the SW differential that exhibits poles at the branch points. This 
contribution can be eliminated at the cost of introducing e-corrections to the SW relations. 

So far a rigorous derivation of the action cut-off term has not been provided. As for the ^-ensemble, it could arise 
from a study of loop equations. These specific Schwinger-Dyson equations are obtained by exploiting the invariance 
of the measure of integration. Those equations are known to be related to the invariance of the Liouville correlator, or 
the /3-ensemble partition function, under the generators of a Virasoro algebra. Similar equations can be established 
for the Nekrasov partition function from its expression as contour integrals. Moreover, Virasoro generators seems 
to act on the Young tableaux summations by adding or subtracting boxes, as suggests the study performed at P = 1 
in [53]. It would be interesting to understand how this interpretation may show up in our formalism. Another major 
motivation to develop the loop equation technique for Nekrasov partition functions is the possibility to access to £2 
corrections. Like for the matrix models, it must be possible to derive a tower of loop equations that can be solved 
recursively, possibly employing a generaUzation of the Eynard-Orantin topological recursion. We hope to address 
these issues in a near future. 

In the second section, we presented the AGT relation as a change of variable between quantum curves. Starting 
from the Schrodinger equation that characterized the Liouville four points function in the semiclassical fimit, we 
derive the Baxter relation proposed in [27] and studied in the first section of the paper. The transformation between 
wave functions was deduced and analyzed at first orders in e. Consistency with the definition of SW differentials 
was also observed. 

The study presented in this paper is restricted to the case of a single SU{Nc) gauge group. It would be tempting 
to extend the picture to the whole set of quiver theories, as was done in [24] for the SW limit, lately extended to the 
NS background [28]. To each node of the quiver diagram should be associated a different density p, or limit shape 
/. On the other hand, densities introduced in the /^-ensemble representing Toda correlators are in correspondence 
with the number of fields, i.e. the rank AT^ — 1 of a single group. The generalization of the transformation between 
wave functions associated to each density is still an open problem. 
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A Definitions and properties of the functions 7 

The function ■~iti,t2 is defined as the logarithm of the Barnes double gamma function [6]. It is invariant under the 
exchange of the parameters e\ and 62, and can be expressed as the following integral, 



7ei,e2(x) = eie2logr2(a;|ei,e2) = 6162 ^ 



1 dt ^ e-*^ 

-rt^T, — 7T7^ — 7T- (A.1) 



We included here the factors €162 to have well defined limits 



€2-^0 as 



1 f°° dt e~*^ 
W)L T^'^Jl^y 7(x) = lim7.(x) = --(logx-3/2). (A.2) 



s=0 

The 7-functions satisfy the difference equations 

(1 - e'^^-)^,,,,,{x) = -e27;(x), (1 - e^^-)7.(x) = -e^'{x). (A.3) 
The expansion at large x for | arg(x)| < tt of the function ^ei,e2 

{x) is [54] 

7ei,.2 W = -\x\\ogx - 3/2) + ^(ei + e2)x(logx - 1) - ^ (^ei€2 + \{A + logx + 0(1) (A.4) 
which implies 

7,(x) = -^x2(logx - 3/2) + lex(logx - 1) - ^e' logx + 0(1). (A.5) 

This expansion is consistent with the limit e ^ reproducing 7(0;). The integral expression for 7^ is not singular 
at t = and can be computed without the (^2-regularization. It is proportional to the derivative of the digamma 
function with argument xje. Integrating twice, and comparing the asymptotic expansion of 7^ with log F, we get 

7^(x) = -elogr[x/€] -xloge + |log(27re). (A.6) 

Instead of 7^, it reveals simpler to study directly the function 7^ defined through its relation to the function G^^ 
given in (2.29), 

eG^\x) = 2(1 - e-'^'^){l - e'^^)^,{x). (A.7) 

This definition allows for an arbitrary shift of 7^ by an affine function of x. Exploiting the relation (2.34) between 

G^^ and 7e, we choose 

la e 
%{x) = + he{x) - -xlog2 + cste, (A.8) 

where the constant is left undefined as it plays no role in our discussion. ExpUcitly, 

7^(x) = -elogr[a;/e] - (e/2) log x - xlog e + (e/2) log(7re). (A.9) 

From the property 

r[dr[-xl = — , (A.10) 

X sin TTX 

we derive the relation 

7g(— x) = — 7e(x) + el{x/e), with l{x) = log | sin7rx|, (A.ll) 

up to shifts of ±Z7re/2 and itiTre, depending on the position of x with respect to the branch cut of the logarithm, and 
on the sign of the sine. 
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B Analysis of a simple difference equation 

In this appendix, we study the difference equation 

(l-e-'^*^ p{t) = S{t-a). (B.l) 

This equation admits solutions defined up to a translation by an e-periodic function. It has a formal solution in the 
form of an infinite series of S functions. However, we instead perform here a perturbative analysis in e. At the first 
order, we find 

ep'oit) = S{t - a) po{t) = ^e{t - a) + Jo (B.2) 

where @{t) is the Heaviside function, and 70 a constant of integration. This constant can be fixed by requiring that 
po has a finite norm using the principal value regularization. We find the value 70 = — l/2e which gives 

poit) = ^sign {t - a), [ po{t)dt = --. (B.3) 

At the next order, decomposing p = po + epi, we find 

Piit) = ^J{t-a). (B.4) 

We choose here the integration constant to be zero in order to have a finite norm for p. We shall take the same choice 
for all higher orders. Up to the fifth order, we have 

p{t) = Isign {t-a) + ^6{t -a) + ^e5'{t - a) - l^e^S^^\t - a) + 0{e'). (B.5) 

At the order O(e'), the density pi is given by the sum over the {k — l — l)th derivative of pk, which are proportional 
to S^^~^^ {t — a) by recursion. It impUes that the density p is expressed as the infinite series 

p{t) = -sign {t-a) + Y^ e^-'ci5^^-'Ht - a). (B.6) 
^ 1=1 

The coefficients q can be obtained from the small x expansion of the generating function x/{l — e~^). We find 
Co = 1, ci = 1/2 and for / > 1, q is vanishing when / is odd, and is expressed in terms of the Bemouilli numbers 
q = Bi/U when I is even. It follows from the expression (B.6) that the moments 



Mfe= / t''p{t)dt (B.l) 
J -A 

can be computed using only the k + 2 first order terms po, - ' ' i Pk+i in the e-expansion of p. In the definition of the 
moments /Hfc, we used a principal value regularization, and A should be sent to infinity at the end of the calculations. 
Using this perturbative approach, we find the following values of the moments, 

A2 1 

epo = -a + e/2, ^1^^ ^ ~2 ~ 2 + ^^/^) > 

Cut-off contributions only appear for even moments. 
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There exists a different approach to compute the moments based on the following formula, 

p{t) (l - e'^A t^dt = a^- ^ \{e - + (e + A)'=+^ - A*+^ - {-Kf+^ 



(B.9) 



This expression is established using a change of variable over t to transpose the shift operator over the density p in 
order introduce the difference equation (B.l). The second term is a boundary term coming from the finite bound of 
integration used in the regularization. In the process, the delta functions (and derivatives) evaluated at itA = a and 
zbA = a — e are discarded since A is much larger than a. Note that the boundary term is independent of a, and will 
be canceled when considering the difference between ppen. and p„, . Such a boundary term typically occurs when the 
density p is replaced by the shape function /. At first orders, we have 

j p{t) (l - e^^*) = 0, j p{t) (l - e^^*) tdt = a- e/2, J p{t) (l - e^^*) t^dt = - - 

(B.IO) 

In particular, one can check the compatibility between the two methods, 

j p{t) (l - e^^*^ tdt = -eno, j p{t) (l - e'^^ {f - et)dt = -2epi. (B.ll) 

To make contact with the work described in [37], we now employ a Fourier transform to solve (B.l). In the 
Fourier space, 

p(t) = / pit)e^'^dt = - ] ,^, e'-(°-^/^). (B.ll) 
Jffi 2ism(eT/2) 

By linearity, it is sufficient to take the sum over the color index to solve the equation (2.24) for Ppert.- The expression 
for the Fourier transformed of the primitive ppert. is obtained after dividing ppert. by it, and writes 

1 ^= 

We recover here the expression of pg^ = — Ppen. found in [37] upon multiphcation above and below by cos(eT/2). 
The pole at r = is related to the norm of ppen. (t) which is divergent. However, this infinity is compensated by the 
presence of p„j in p and we do not see any reason here to impose a finite norm, in contrast with the claim of [37]. 

C Derivation of the change of variable in the massive case 

Starting from the Schrodinger equation (3.12), we again multiply on the left by x{x — l){x — q). It ehminates the 

double poles at .t = 0, 1, leaving only single poles. The remaining poles at a: = 1 and x = qwe eliminated by a 
good choice of the function s{x) in the change of variable z = xy + s{x). With this change of variable, we have 

a;2y2 = ^2 _ 2zs{x) + Exs' {x) + s{xf - e{z - s{x)). (C.l) 
Using commutation relations, we push all x-dependence to the right, and write the equation as 

(^x Z{x) + R{x)x-^ -{q- l)E = 0, (C.2) 

In this expression, terms that may present poles at x = 1 and x = q have been gathered in the operator R{x) which 
is ^-independent, 

Rix) = 2es(x)(x2 -q) + (x-l){x- q) [sixf + exs'ix)] - , ^^'^""^ -. (C.3) 

[x — l){x — q) 
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The remaining z-dependent term is free of singularities at x = l,q provided that s{x) has no more than first order 
poles, 

Z{x) = {z- e){z - 2e)x - (1 + q)z{z - e) + qz{z + e)x~^ + (e - 2z)s{x)x~'^{x - l){x - q). (C.4) 

The function s{x) is determined from the requirement that R{x) is a polynomial in x of degree at most two. In 
particular, poles atx = l,q must vanish. We propose the ansatz 

s{x) = a + ^^ + ^^. (C.5) 
x — 1 X — q 

To fulfill the requirement of vanishing residue at x = 1 and x = q on R{x) we take si = —mi and Sq = —qm2. 
It is worth noticing that the constants s\ and Sq are solutions of a quadratic equation and should a priori be infinite 
series in e. Surprisingly, the solution we find has no e-corrections. Under this choice of s{x), and regrouping the 
terms according to their a;-dependence, we find 

f^^j^x - (1 + q)P{z) + r{z + e)x^^ 
r{z) 



^ = 0, (C.6) 



which is exactly the form of the Baxter TQ relation (2.66). In general r{z) depends on three variables mo, mi and 
1712, which leads to an incorrect form of the function M{z). To reproduce the correct mass polynomial as given in 
(3.1), we have to set a = —mi — e/2 in s{x) which suppress the dependence in the mass mi of r{z). Then, the 
zeros of r{z) coincide with the values of /x+ and /x_ given in (3.5), and we recover the expression (3.1) for M{z). 
Finally, the polynomial P{z) is indeed monic and of degree two,^^ 

2q 

= - (ei + 62)2; + 6162, ei + 62 = — — (m2 - m-i), 

J- "r 9 (C.8) 

(1 + g)eie2 = —(1 — q)E — q (m^ — 3m| + 2(mi + e)m2 — m\ — mg) — 2m2(m2 — e) — 2mQ. 

Gathering all the constraints, the function s{x) can be expressed with the help of the /3-ensemble potential as in 
(3.21). In the massless Umit, it reduces to s{x) = —e/2. 



D Study of the Mellin transform between wave functions 

In this appendix, we examine the quantum change of variable z = xy — e/2 for the wave functions ipa,'y = 
2,a/eg-7x/e supposc here 7 > 0, in the case 7 < the integration range [0, +00 [ of the Mellin transformed 
should be replaced by ] — oo,0]. These wave functions satisfy a Schrodinger equation with a potential having a 
double pole at x = 0, 

y'^^Pa^ix) = T{x)^a'yix), T{x) = a{a - e)^ - 2a^^ + 7^. (D.l) 

Such a potential would appear in the study of a sphere with a single puncture at x = 0. Performing the quantum 
change of variable over this Schrodinger equation with r{z) = z^ — [a — e/2)^, we get the difference equation 

-i^M{z)i^o.^{z - e) - 2a7V'a7(z) - V^^^l^ + e) = 0, M{z) = r{z). (D.2) 

'''Another solution is possible, s\ — mi — e and Sq = qijUq — e). In this case, to identify properly the mass polynomial M{z), we need to 
take into account the action of the Weyl group of the SU (2) gauge symmetry that flips a — >■ —a, mo — >■ —mo, moo — >■ —moo, rtii t — m\ 
and m2 — >• e — m2 [6]. It leads to define 

/i± = e — m2 ± mo, /i± = e — mi ± moo, (C.7) 
and we recover the correct Baxter relation with a different polynomial P{z). 
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It is easy to work out the transformation of tpa-y, we get 

^ (z) = {z-a)/e+l/2^{z+a)/e-l/2-^^y2 + {z + a)/e]. (D.3) 

cos((2; — aj/e) 

It is also trivial to check that this wave function indeed satisfies the difference equation (D.2). 

It was shown in [35] that the definition of the SW differential coincides with the coordinate transformation that 
trivializes the stress energy tensor of the Coulomb gas representation. In these new coordinates, that for simplicity 
we still denote x and y, we have dS = dx and T{x) = 1. This problem corresponds to the case we just treated, 
with the special values a = and 7 = 1. We note that the difference equation simplifies in this case since the term 
proportional to ip{z) disappear. Under the choice of coordinates r{z) = z + e/2, the factor M{z) disappear and the 
trivial wave function = 1 solves the equation. This choice of coordinate modifies the integration kernel of the 
Mellin transform through the function a{z) and indeed produces a trivial, i.e. e-periodic, wave function. 
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